########## MAURCIO TAKASHI COUTINHO WATANABE ########### Exerccio 8

getwd()
setwd("C:/Users/Maurcio/Documents/aula8")


##### Galileu estava certo?
# Utilizando os dados disponveis no tutorial para verificao do melhor modelo (se o de 2 ou 3 grau)
init.h = c(600, 700, 800, 950, 1100, 1300, 1500)
h.d = c(253, 337, 395, 451, 495, 534, 573)
mod2 <- lm(h.d~init.h+I(init.h^2))
mod3 <- lm(h.d~init.h+I(init.h^2)+I(init.h^3))
anova(mod2,mod3)
# Um polinmio de 3 grau parece ser um modelo melhor do que o do de 2 grau tendo em vista que ele apresenta um menor valor de RSS.



##### Massa de Recm-Nascidos
babies <- read.table("babies.txt",header=TRUE,sep="")
babies

babies <- babies[babies$bwt!="999",]
babies <- babies[babies$gestation!="999",]
babies <- babies[babies$parity!="9",]
babies <- babies[babies$height!="99",]
babies <- babies[babies$weight!="999",]
babies <- babies[babies$smoke!="9",]
babies <- babies[babies$age!="99",]
babies

# No foram comparadas as variveis parity e smoke (no h muito sentido em analis-las, j que so praticamente "binrias")
mod.a <- lm(bwt~gestation, data=babies)
mod.b <- lm(bwt~age, data=babies)
mod.c <- lm(bwt~height, data=babies)
mod.d <- lm(bwt~weight, data=babies)
anova(mod.a,mod.b,mod.c,mod.d)
# o melhor modelo  aquele relaciona massa ao nascer e o tempo de gestao

# Agora vamos adicionar o efeito de outras variveis
mod.e <- lm(bwt~gestation+age, data=babies)
mod.f <- lm(bwt~gestation+height, data=babies)
mod.g <- lm(bwt~gestation+weight, data=babies)
anova(mod.e,mod.f,mod.g)
# Podemos perceber que a adio da varivel altura da me, melhora o modelo.

# Vamos agora verificar se encontramos um modelo melhor adicionando as outras variveis
mod.h <- lm(bwt~gestation+height+age, data=babies)
mod.i <- lm(bwt~gestation+height+weight, data=babies)
anova(mod.h,mod.i)

mod.j <- lm(bwt~gestation+height+weight+age, data=babies)
anova(mod.i,mod.j)

# Podemos perceber que mesmo adicionando outras variveis, no h um ganho muito significativo a ponto de justificar suas incluses.
# Neste caso o modelo que relaciona a massa dos recm nascidos ao tempo de gestao e a altura da me (mod.f) apresenta o melhor "custo-benefcio".